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The electronic band structure of GaTe has been calculated by numerical atomic orbitals density- 
functional theory, in the local density approximation. In addition, the valence-band dispersion along 
various directions of the GaTe Brillouin zone has been determined experimentally by angle-resolved 
photoelectron spectroscopy. Along these directions, the calculated valence-band structure is in 
good concordance with the valence-band dispersion obtained by these measurements. It has been 
established that GaTe is a direct-gap semiconductor with the band gap located at the Z point, that is, 
at Brillouin zone border in the direction perpendicular to the layers. The valence-band maximum 
shows a marked p-like behavior, with a pronounced anion contribution. The conduction band 
minimum arises from states with a comparable s- p-cation and p-anion orbital contribution. Spin- 
orbit interaction appears to specially alter dispersion and binding energy of states of the topmost 
valence bands lying at F. By spin-orbit, it is favored hybridization of the topmost p^-valence band 
with deeper and flatter Px-Py bands and the valence-band minimum at F is raised towards the Fermi 
level since it appears to be determined by the shifted up Px-Py bands. 



FAGS numbers: 71.20.Nr, 79.60.Bm 
I. INTRODUCTION 

Layered materials of the III- VI family like GaSe, InSe, 
and GaTe are of special interest for their potential .-ap- 
plications in nonlinear and optical bistable devicesJil'B as 
well as iiVilAe development of solar cells and solid-state 
batteries. ou Among these compounds, GaTe is one of 
the less studied because it presents a higher anisotropic 
and more complex crystal structure than the rest of the 
III- VI materials. GaTe belongs to the B2/m space group 
and has a base centered monoclinic unit celLQ The layer 
structure [see Fig. ||(a)], as that of GaSe or InSe, is com- 
posed of a four-sheet anion-cation-cation-anion intralayer 
stacking pattern in which the bonds are mainly covalent 
with some ionic contribution. Nevertheless, in contrast 
to GaSe and InSe, one-third of the cation-cation bonds 
in GaTe lies almost in the layer plane in a direction per- 
pendicular to the c axis. Differences also exist when con- 
sidering layer stacking sequences. Similarly to the rest of 
III- VI compounds, interlayer anion-anion bonds in GaTe 
are of weak van der Waals type. Generally, this fact 
promotes the existence of polytypes with different-stack- 
ing sequences, as is observed in InSe and GaScffl Nev- 
ertheless, no polytypism seems to occur in GaTe, as is 
suggested by the absence of conjugate modes in Raman- 



scattering and infrared-absorption measurements J13 The 
whole structure gives rise to fourfold coordination for the 
Ga atoms, with three Te and one Ga atoms, and to three- 
fold coordination for the Te with Ga atoms. As a con- 
sequence of these differences in the crystal structure of 
GaTe with respect to other layered compounds, GaTe 
has only one twofold symmetry axis, lying along the c 
axis, whereas other III- VI compounds have a threefold 
symmetry axis perpendicular to the layer plane. 

Besides the common in-plane-out-of-plane anisotropy 
observed in the optical and electrical properties of the 
layered compounds, the special in-plane anisotropy in 
bonding in GaTe makes the optical properties also 
anisotropic in the layer plane. GaTe is a direct-gap semi- 
conductor wi|th strong excitonic absorption even at room 
temperatureJlJilJ The optical properties near the banck 
gap energy region appear to be slightly anisotropic.E£l 
Anisotropy of the absorption coefficient has been ob- 
served to increase at higher energies, which was at- 
tributed to transitioaa.from deep states related to the in- 
plane Ga-Ga bondso Moreover, low-temperature pho- 
toluminescence and transmission measurements revealed 
the presence of different excitpa states depending on 
the in-plane polarization vector .Ej Recently, these results 
were explained by the existence of three types of exciton 
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states due to j-j coupling, whose optical transitionp-are 
allowed under the group symmetry of this materialJlj 

There is still a lack of both theoretical and experi- 
mental studies on the electronic band structure of GaTe. 
One of the reasons for this lacuna lies on the complex- 
ity of the primitive unit cell, that contains 6 GaTe units 
(54 valence electrons), which has made band structure 
calculations hardly available up to now. Recently, the 
electronic band structure of GaTe has been calculated 
along some directions by the_a& initio tight-binding linear 
mufhn-tin orbitals method.cj These results are in agree- 
ment with the direct-gap behavior of this semiconductor 
and locate the band edge at the border of the Brillouin 
zone (BZ). Nevertheless, there are still some open ques- 
tions regarding the electronic band structure along the 
high-symmetry directions of the GaTe BZ, as well as the 
experimental confirmation of the electronic band struc- 
ture obtained in the calculations. 

In this work we study the electronic properties of the 
GaTe compound. The GaTe band structure has been cal- 
culated by numerical atomic orbitals density-functional 
theory (NAO-DFT). In addition, angle-resolved photoe- 
mission (ARPES) measurements have been carried out 
in GaTe(OOl) single crystals. From these measurements, 
the valence-band dispersion of GaTe was extracted and 
interpreted on the basis of the band structure calcula- 
tions. The experimental setup is described in Sec. II. 
Section III is devoted to the band structure of GaTe cal- 
culated along different directions of the GaTe BZ, includ- 
ing the technical aspects of the calculation. In Sec. IV, 
the experimental results obtained by ARPES are pre- 
sented and discussed in the frame of the calculated GaTe 
band structure. Finally, in Sec. V we present the main 
conclusions of this work. 



II. EXPERIMENTAL 

The experiments were performed at LURE (Orsay, 
France) using the Spanish-French (PES2) experimen- 
tal station-. of the Super- Aco storage ring, described 
elsewhere. 113 The measurements were carried out in a 
purpose-built ultra-high vacuum system, with a base 
pressure of 5xl0~^^ mbar, equipped with an angle re- 
solving 50 mm hemispherical VSW analyzer coupled on 
a goniometer inside the chamber. The manipulator was 
mounted in a two-axes goniometer which allows rotation 
of the sample in the whole 360° azimuthal angle and in 
the 180° polar emission angle relative to surface normal 
{Qoff), with an overall angular resolution of 0.5°. Pho- 
toelectrons were excited with p-polarized synchrotron ra- 
diation in the 20-60 eV energy range. For this selected 
photon energy (hu) range, the energy resolution was of 
60 nieV. With this setup, different energy-distribution- 
curve (EDC) series were recorded scanning the hz/ as well 
as 6o//, at constant incident angle of the light (6i). In 
these conditions, changes of polarization effects on ini- 



tial states along one series are neglected. All the EDC 
spectra shown in this work are referred to the Fermi level 
(Ei.). 

The GaTe single crystals used in this work were 
cut from an ingot grown by the Bridgman-Stockbarger 
method. The samples were cleaved in situ after intro- 
duction in the ultra-high vacuum chamber. The samples 
were easily cleaved in the layer plane due to the existence 
of weak inter layer van der Waals bonds. LEED spectra 
showed a sharp spot pattern corresponding to the bulk 
monoclinic material. No surface impurities were detected 
by synchrotron radiation photoemission measurements. 
The samples were oriented by azimuthal and polar photo- 
electron diffraction scans recording the Ga 3 c? peak inten- 
sity. Unambiguous azimuthal orientation of the samples 
was possible due to the B2/m symmetry of the crystal. 



III. BAND STRUCTURE OF GATE 

A. Calculation details 

In this work, w& present fully self-consistent density- 
functional theoryli3 (DFT) calculations of the electronic 
structure of GaTe. The calculation^ are performed 
in the local-density approximationEj The exchaiige- 
correlation potential is that of Ceperlejt-and AdlerE^I as 
parameterized by Perdew and Zunger.tj Only the va- 
lence electrons are considered in the calculation, with 



the core being replaced 
ativistic pseudopptentialf 



norm- conserving scalar rel- 
factorized in the Kleinman- 
Bylander form.tJ The pseudopotentials were generated 
using the atomic valence configurations 5s^5p^ for Te 
and 3rfi°4s24pi for Ga. The cutoflf radii were 2.0, 2.0, 
3.0, and 3.0 a.u. for the s, p, d, and / components in Te, 
respectively, and 2.1, 2.3, 1.1 and 2.0 a.u. for the s, p, 
d, and / components of Ga, respectively. Since the core 
and valence charges overlap significantly for both atornic 
species, we include nonlinear partial-core corrections^ 
with matching radii of 1.0 a.u. for Te and 0.5 a.u. for 
Ga to describe the exchange and correlation in the core 
region. 

The valence one-particle problem was solved using a 
linear combinatiotupf numerical (pseudo-) atomic orbitals 
with finite range.E3 The details of the basis generation 
(including multijple-C and polarization functions) can be 
found elsewhere Hj In this work we have used a split- 
valence double-^ basis set with a single shell of polar- 
ization orbitals (that is, containing two s shells, two p 
shells, and one d shell both for Tc and for Ga), as ob- 
tained "adth an energy shift of 250 meV and a split norm 
of 15%.l3 The integrals of the self-consistent terms of the 
Kohn-Sham Hamiltonian are obtained with the help of a 
regular real space grid on which the electron density is 
projected. The Hartree potential is calculated by means 
of fast Fourier transforms in that grid. The grid spac- 
ing is determined by the maximum kinetic energy of the 
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plane waves that can be represented in that grid. In the 
present work, we used a cutoff of around 300 Ry (which 
changes slightly with the volume of the unit cell) . A reg- 
ular grid of 5x5x5 k pointsEZI was used to sample the BZ. 
We have checked that the results are well converged with 
respect to the real space grid, the BZ sampling and the 
range of the numerical atomic orbitals. _The calculations 
were performed using the Siesta code.t3 

In our band structure calculations, we used a centered 
monoclinic cell with primitive vectors a' = — b-|-(c— a)/2, 
b' — — b — {a + c)/2, and c' — — b, where a, b, and c are 
the repeat vectors of the monoclinic cell [see Figure |l|(a)] . 
Note that, with this election, a' and b' vectors are in the 
layer plane. 



B. Brillouin zone and electronic band structure 

Figure |l|(a) shows the crystal structure of GaTe, de- 
scribed previously. In this figure, we have indicated the 
crystallographic axes of the GaTe monoclinic unit cell. 
The GaTe crystal structure is quite similar to those of 
GaSe and InSe, except for the existence of cation-cation 
bonds lying almost in the layer plane, with the same co- 
ordination and bonding behavior. The comparison can 
be also extended to the reciprocal space. Figure |l](b) 
shows the calculated GaTe BZ, in which the main points 
are indicated. The relation between the different symme- 
try directions of the GaTe BZ and the reciprocal lattice 
vectors of the centered monoclinic unit cell is shown in 
Table |. 

The GaTe BZ can be easily described in terms of the 
crystal structure. The binary axis corresponds to the FY 
high-symmetry direction and a crystal vector perpendic- 
ular to the layer plane defines a direction parallel to the 
rZ high-symmetry direction, in the reciprocal space. The 
GaTe BZ resembles that of InSe, but contracted along 
the FX direction. The contraction is due to the pres- 
ence of the quasi-in-plane Ga-Ga bonds in GaTe, which 
produces that the fundamental translation in the layer 
plane consists of three Ga2Te2 unities, in contrast to one 
unity in the other III- VI compounds. On the other side, 
as both GaTe and InSe have only one layer per unit cell, 
one should expect the electronic band structure of GaTe 
to be related to that of InSe, after a three-folding in the 
layer plane. 

Figure || shows the electronic band structure of GaTe 
along different directions of the GaTe BZ [Fig. |^(b)] cal- 
culated by the NAO-DFT method described previously. 
The origin of binding energy has been taken at the Ej;^. 
In spite of the large number of bands obtained, a de- 
tailed description of the valence band can be done. The 
GaTe unit cell consists of three Ga2Te2 unities. There- 
fore, there are per unit cell: three Ga-Ga bonds, 18 Ga-Te 
bonds, and one electronic pair associated to each one of 
the six Te. This gives rise to 27 doubly occupied bands, 
which are precisely those appearing with binding energies 



between 14 eV and E p and composing the GaTe valence 
band. Below this energy (at 16.5 eV), the 3d-core lev- 
els of Ga are present (not shown). The unoccupied side 
of the band structure is rather more complicated to be 
described, since, besides the 18 Ga-Te and 3 Ga-Ga an- 
tibonding levels, all contributions from the polarization 
d-orbitals overlap. 

A detailed band to band description of the band struc- 
ture of GaTe would not be very useful due to the fact 
that the low symmetry of the crystal makes the orbital 
mixture of each band to be very high. Nevertheless, a 
general description of the orbital character of the GaTe 
band structure can be performed by means of the density 
of states (DOS). Figure ||(a) shows the total DOS calcu- 
lated by the NAO-DFT method. Figure |(b) and |(c) 
show the different projections of the s, p, and d orbitals 
onto the DOS for Ga and Te, respectively. These re- 
sults indicate that the valence band of GaTe can be di- 
vided into three different groups: (i) One group of six 
bands appearing between 14 and 11 eV, which show a 
pronounced Te 5s relative behavior with a small Ga 4s 
and Ap contribution, (ii) a second group of six bands ap- 
pearing between 8 and 5 eV, which mostly show a Ga 
4s character with a pronounced Te 5p contribution in 
the low binding energy side, and (iii) a third group of 
15 bands comprising the range of energy between 5 eV 
and the E^?, which are a mixture of Te bp and Ga 4p 
orbitals with a small Ga 4d and 4s contribution. The 
origin of these groups differs from one to the other. On 
one side, the first and second group of bands appear to 
be mostly originated by the Ga-Te bonds, showing those 
with a marked p-like behavior a more pronounced dis- 
persion along the FZ high-symmetry direction (as well as 
along directions perpendicular to it) than those with a 
predominant s-like behavior (Fig. On the other side, 
in the third group of bands, they coexist six bands orig- 
inated also from the Ga-Te bonds, three bands from the 
Ga-Ga bonds, and the six electronic pairs associated to 
the Te (mostly Te 5p). Focusing our attention on the 
topmost valence bands, the fact that their orbital char- 
acter has an important p-like component perpendicular 
to the layers, makes that most of the bands display a 
noticeable dispersion along FZ high-symmetry direction 
(Fig. I). 

With regard to the DOS corresponding to the unoc- 
cupied side of the band structure, the huge overlap of 
Ga-Te and Ga-Ga antibonding levels with the d-orbitals 
of polarization makes an orbital description quite more 
complicated than that carried out for the occupied side 
of the DOS. Nevertheless, it should be mentioned that 
states at the conduction band edge are mostly originated 
from Te bp and Ga Ap orbitals, with an important con- 
tribution from Ga 4s and a smaller contribution Te 5d 
and Ga Ad, which suggests a pronounced spz-like behav- 
ior to the conduction band edge. Only for higher states 
in the conduction band, contributions from Ga and Te 
rf-orbitals of polarization appear to be of the same mag- 
nitude than those of the p-orbitals. 
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The general shape and the different orbital contribu- 
tions to the calculated DOS by NAO-DFT are quite sim- 
ilar to those obtained previously by the ab initio tight 
binding mcthodEj Some differences can be appreciated 
in the unoccupied side of the DOS. By both methods, the 
behavior of the conduction band edge seems to be estab- 
lished as a mixture of s and p orhitals. Nevertheless, in 
contrast to previous calculationsjlj the Ga 4p contribu- 
tion to the conduction band edge appears to be of the 
same order as that of the Ga 4s, which gives a more 
pronounced p-like behavior to the conduction band edge 
in the present calculations. Also the contribution of Ga 
and Te d-orbitals of polarization to the conduction band 
edge, as well as to deeper states of the conduction band, 
seem to be underestimated in those calculations (mostly 
the Ga Ad). As regards the band dispersion, authcpj 
seem to have used a standard six-face monoclinic BZ,E£l 
which does not correspond to the GaTe BZ. The band dis- 
persion calculated along some directions in Ref. ^ (see 
Table |) is very similar to the present results. However, 
along other directions, as the FZ and AM directions in 
Fig. 8 of Ref. the path crosses out of the first BZ 
to neighboring ones, giving rise to the gitasi-symmetrical 
band dispersion observed along these directions. 

From the above results, it can be concluded that GaTe 
is a direct-gap semiconductor with the gap located at 
zone border at Z point of the GaTe BZ, which has the 
same symmetry as F (B2/m). The band gap appears to 
be underestimated (1.098 eV), which is a well-known ten- 
dency o£-the LDA approximation in semiconductors. As 
in InSe,E3 the valence-band maximum is mostly formed 
by antibonding n^orbitals of the anion. Nevertheless, in 
contrast to InSeJ^ the conduction-band minimum is not 
only formed by the s-orbital contribution of the cation. 
In GaTe a more intense contribution of cation and anion 
p-orbitals can be observed. This can be attributed to the 
fact that the presence of the quasi-in-plane cation-cation 
bonds distorts the anion disposition and, therefore, the 
interaction of Te atoms with those of contiguous layers 
and with Ga atoms increases. 



IV. EXPERIMENTAL RESULTS AND 
DISCUSSION 



In this section, we compare the calculated electronic 
band structure of GaTe with the band dispersion de- 
termined experimentally by ARPES in GaTe(OOl) single 
crystals. First, we will show and discuss the electronic 
band dispersion along the FZ high-symmetry direction, 
proceeding later with those determined along different 
directions perpendicular to it. 

Figure ^(a) shows an EDC series recorded at normal 
emission with constant 8^ in a GaTe(OOl) sample, as a 
function of the hv. In these spectra, several peaks have 
been identified and labelled by small solid bars. In some 
of them, two sharp dispersing peaks can be also observed. 



which can be attributed to second-order transitions from 
deeper occupied states induced by the experimental sys- 
tem. As the hi/ increases, these sharp peaks move to 
lower binding energies. Therefore, photoemission sig- 
nal from deeper valence-band transitions than those of 
Fig. U(a) can be studied in EDCs recorded by extending 
up the hj/-range of Fig. ^(a). Figure, ^(b) shows the re- 
sults of these measurements. In these spectra, only broad 
features can be identified, at binding energies of 12.6 eV, 
which have been labelled as s-band. 

By these results, the dispersion of states along the FZ 
high-symmetry direction can be determined. Neverthe- 
less, to compare these results with the calculated band 
structure, the momentum of photoelectrons should be 
estimated. Photoemission process involves energy and 
crystal momentum conservation. In the frame of the 
three-step model of the photoexcitation mechanism in 
ARPES, both momentum components, and of the 
photoelectrons inside a bulk material can be expressed 
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BE~<^)cos^eoff)-Vo , (lb) 



where BE is the binding energy, $ is the work function, 
Vo is the inner potential, m is the free-electron mass, 
and h is the reduced Planck constant. A complete de- 
termination of the photoelectron momentum requires an 
accurate estimate of the Vo [Eq.(lb)]. We have assumed 
it as a constant parameter. An initial testing value of 
Vo^-l9 eV ($~6 eV) can be obtained from the bottom 
of the valence band (Fig. ||). After testing several values 
of Vo and taking into account that krz=0.421 A~^, we 
have adopted a value of Vo=-25.'LeV, which is two times 
higher than that found for InSe.Ell Such a high value of 
the Vo for GaTe is not surprising, since the orbital hy- 
bridization giving rise to the GaTe valence-band struc- 
ture is higher than in InSe. In any case, it is reasonably 
close to the value expected from the valence-band calcu- 
lations. 

Figure ^(c) shows the valence-band dispersion ex- 
tracted from the peaks identified in Fig. |(a), by means 
ofEqs. (U). The calculated valence-band dispersion along 
the FZ high-symmetry direction has been also included. 
In spite of the large number of different bands present, 
it can be observed that the calculated band structure re- 
produces most of the traces observed in the experimen- 
tal band dispersion. The large concentration of quite 
flat bands appearing at lower binding energies is well re- 
produced by the calculated valence-band structure, with 
states lying at similar binding energies. These facts also 
put forward the highly hybridized behavior of top states 
of the valence-band structure predicted by calculations 
(Fig. H). Besides this, dispersion of states with a marked 
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Te 5p-Ga As orbital character (that is, states between 8 
and 5 eV in binding energy) is quite well reproduced. 
Moreover, the presence of the broad s-hand peaks in 
Fig. ^(b) can be attributed to transitions from the Te 5s- 
like states, located at binding energies of ^12 eV (Fig. |^). 
Nevertheless, some discrepancies appear for states ly- 
ing at the top of the valence band. Experimentally, the 
valence-band maximum along the FZ high-symmetry di- 
rection appears to be at Z point and the minimum of the 
valence-band lies at F. These facts are in concordance 
with the calculated valence-band structure. Moreover, 
dispersion of the topmost band around Z is well repro- 
duced by the calculations. However, the valence-band 
minimum at F obtained by calculations appears shifted 
down by 0.8 eV with respect to the experimental one. A 
similar situation occurs with the second topmost valence 
band. For this band, the calculated band minimum, also 
at F, appears to be shifted down by 0.6 eV with respect 
to the experimental one. This disagreement will be dis- 
cussed later on. 

Figure ||(a) shows an EDC series measured by ARPES 
with hv=?>Z eV in a GaTe(OOl) sample. The sample ori- 
entation has been chosen to probe, by these ARPES mea- 
surements, points of the reciprocal space contained in the 
FZH plane. In fact, for this selected hv, these points are 
expected to have a kj^ lying quite close to the Z point 
[Ecs. (^] . In the spectra shown in Fig. ||(a) , several peaks 
have been identified and labelled by small solid bars. Fig- 
ure |^(b) shows the band dispersion extracted from the 
peaks identified in Fig. ||(a). 

In order to compare the experimental band dispersion 
shown in Fig. ||(b) with the calculated band structure, it 
should be taken into account that, by the ARPES mea- 
surements summarized in Fig. ^(a), particular points of 
the reciprocal space along a parabolic curve are probed 
for each constant initial state [Eqs. (|l|)]. For initial states 
with low binding energy, the minimum of the parabola 
will be close to the Z point, but when running from Z to 
H (k^H =0.744 A^^) the parabolic curve perpendicularly 
deviates from horizontality by 0.07 A~^. For deeper ini- 
tial states, the parabolas probe points somewhere along 
the FZ high-symmetry direction. These facts produce 
that, by these measurements, states along the ZH di- 
rection are not strictly probed. Nevertheless, it has been 
found that dispersion of valence-band states along the FZ 
high-symmetry direction can be considered negligible in 
a large k_L-range [Fig. ^(c)]. These facts encourages us to 
assume that the calculated band structure along the ZH 
direction is a good approximation for the experimental 
band dispersion, as shown in Fig. ^(b). In this figure, it 
has been also introduced the calculated band dispersion 
of states with binding energies higher than 3 eV along 
the FY high-symmetry direction, which is not far from 
that calculated along the ZH direction and gives a bet- 
ter agreement with the experimental band dispersion of 
deeper bands. 

Figure ^(a) and ^(b) show two different EDC series 
measured by ARPES with hi^=33 cV in a GaTc(OOl) 



sample. In these ARPES measurements, a sample orien- 
tation has been chosen to probe points of the reciprocal 
space contained in the FZM plane. In these spectra, sev- 
eral peaks have been identified and labelled by small solid 
bars. Figure |6|(c) shows the FZM transversal cut of the 
GaTe BZ, in the extended zone scheme. The parabolic 
curve defined by ARPES measurements with hi/=33 eV 
for constant initial states at the Ej? is also included in 
this figure. This scheme illustrates the fact that, by both 
sets of measurements shown in Figs, ^(a) and ^(b), points 
of the reciprocal space are probed along the MZM' and 
MFM' directions, respectively. Figure ^(d) shows the 
band dispersion extracted from the peaks identified in 
Figs, ^(a) and ^(b). Assuming the same hypotheses as 
those discussed previously, we have also included the cal- 
culated band dispersion along the MZM' and MFM' di- 
rections, which reproduce most of the traces of the band 
dispersion obtained experimentally. Focusing again at- 
tention on the two topmost valence-bands, their disper- 
sion around Z along the MZM' direction appears to be 
rather flat, coherently with the behavior expected by 
NAO-DFT calculations. Nevertheless, dispersion of these 
bands along the MFM' direction shows the same behavior 
as that observed along the FZ high-symmetry direction 
[Figs. §(c)], with the same separation at F between the 
calculated and experimental minimums of these bands. 

The above results obtained by ARPES measurements 
along different directions of the GaTe BZ support the re- 
sults obtained by NAO-DFT band structure calculations 
(Fig. H). Among these results, it should be emphasized 
the fact that experimental and calculated band structure 
place the valence-band maximum in GaTe at zone-border 
at Z, similarly as occurs in the InSe compound.Ej Nev- 
ertheless, despite the good agreement obtained through 
the whole reciprocal space studied, there still remains one 
point of discrepancy related to the two topmost valence- 
bands at F [Figs. U(c) and ^(d)], since the binding energy 
of the minimum of these bands appear to be overesti- 
mated by calculations. These facts can be explained since 
spin-orbit interaction was not considered in the present 
calculations. In the related compound InSe, the band 
giving rise to the valence band maximum at Z mostly has 
a marked Pz-orbital character. States of this band have a 
total orbital moment j=l/2. Without considering spin- 
orbit interaction, this band disperses down along the FZ 
direction, showing the minimum at F. Below this band, 
the next ones are quite flat bands with a Px-Py character. 
Spin-orbit interaction splits oS Px-Py states-jwith different 
j, rising those with j=3/2 by 0.3-0.4 eV,E3 with respect 
to the valence band maximum. This produces that the 
minimum of the topmost valence band at F is determined 
by the shifted up Px-Py bands. A more intense spin-orbit 
coupling can be expected in GaTe, since the nature of 
its topmost valence bands is the same as those of InSe 
and its anion has a higher mass. In fact, the spin-orbit 
induced shift on the valence-band minimum at F in GaTe 
may be evaluated. Taking into account that the atomic 
spin-orbit Pzji-Vxji sphtting is of 589 and 247 meV for 
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V. CONCLUSIONS 
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FIG. 1. (a) Schematic view of the crystal structure of GaTe 
and its monoclinic unit cell. Non-equivalent atoms have dif- 
ferent shadowing. Big and small circles correspond to Te and 
Ga atoms, respectively, (b) GaTe centered monoclinic Bril- 
louin zone, in which the main points are indicated. Equiva- 
lent Q-points (Q generic) have been labelled as Q', Q", and 
so forth. The FXY plane is parallel to the layers and the 
FY direction corresponds to the binary axis. The FZX plane 
defines the mirror plane. 



FIG. 2. Band structure of GaTe calculated by the 
NAO-DFT method along various directions of the GaTe BZ. 

FIG. 3. (a) Density of states of GaTe. (b)-(c) Projections 
of the s, p, and d orbitals onto the total density of states, for 
Ga and Te, respectively. 

FIG. 4. (a)-(b) Normal-emission valence-band EDCs mea- 
sured in GaTe(OOl) single crystal, as a function of the hi/. 
The hi/ range scanned has been indicated in both figures. 
The different peaks identified in the EDCs in (a) have been 
marked by small solid bars. The broad peaks identified in 
the EDCs in (b) have been labelled as s-band. (c) (circles) 
Band dispersion diagram extracted from the peaks identified 
in (a). The calculated GaTe band structure along the ZFZ 
high-symmetry direction has been also plotted, in solid lines. 

FIG. 5. (a) Valence-band EDCs measured by ARPES with 
hi/=33 eV in GaTe(OOl) single crystal. For this selected hi/, 
the ZH direction is scanned by ARPES. The different peaks 
identified have been marked by small solid bars, (b) (circles) 
Band dispersion diagram extracted from the peaks identified 
in (a). The calculated GaTe band structure along the ZH 
direction has been also plotted, in solid lines. The calculated 
band dispersion of states with binding energies higher than 
3 eV along the FY high-symmetry direction has been also 
plotted, in dashed lines. As a reference, the fcy-projection of 
some points of the GaTe BZ are indicated at the top. 



FIG. 6. (a) Valence-band EDCs measured by ARPES with 
hi/=33 eV in GaTe(OOl) single crystal. For this selected hv, 
the ZM direction is scanned by ARPES. (b) The same that 
(a) extending the Oo//-range to scan the MFM' direction. 
In both EDC series, the different peaks identified have been 
marked by small solid bars, (c) Schematic view of the BZ cut 
with the FZM plane. The parabolic curve defined by ARPES 
measurements with hi'=33 eV for constant initial states at 
the E_F is also included, (d) (circles, squares) Band dispersion 
diagram extracted from the peaks identified in (a) and (b), 
respectively. The calculated GaTe band structure along the 
MZM' and MFM' directions has been also plotted, in solid 
lines. As a reference, the fcy -projection of some points of the 
GaTe BZ are indicated at the top. 
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